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ABSTRACT: Recently, we found that the correlation between the eigenvalues of ran- 
dom hermitean matrices exhibits universal behavior. Here we study this universal behavior 
and develop a diagrammatic approach which enables us to extend our previous work to 
the case in which the random matrix evolves in time or varies as some external parameters 
vary. We compute the current- current correlation function, discuss various generalizations, 
and compare our work with the work of other authors. We study the distribution of eigen- 
values of Hamiltonians consisting of a sum of a deterministic term and a random term. The 
correlation between the eigenvalues when the deterministic term is varied is calculated. 



I. INTRODUCTION 



We have been studying correlations between energy eigenvalues in random matrix 
theory 1 ' 2 ' 3 ' 4 in an attempt to uncover possible universal behavior in disordered systems. 
Let us begin by summarizing the main results of our previous work.^ Consider an ensemble 
of N by N hermitean matrices ip defined by the probability distribution 

PM = (i.i) 

with 

H(<p) = -^TrV(ip) (1.2) 

for any even polynomial V and with Z fixed by the normalization J dipP(ip) = 1. Define 
the Green's functions 

G <^H^ r ^ tr ^) (L4) 

and so forth, where 

(Ofa)) = J d<pO(<p)P(<p) (1.5) 
The density of eigenvalues is then given by 

P (/jl) = QjtrSfr -<P)} = ^ImG(fi + ie) (1.6) 
and the correlation between eigenvalues, by 

p(w) = (Jj1r6([i-(p)^tr6(v-(p)j 

= (-l/4n 2 )(G(++) + G(— ) - G(+-) - G(-+)) 
with the obvious notation 

G(±,±) = G(ji±ie,v±i6) (1.8) 

(signs uncorrelated) . 

In the large N limit, p{p,,v) — > p{p)p{v), and thus it is customary to define the 
connected correlation 

p c (ji,v) = p(ji,v)- p(jj)p(u) (1.9) 
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Note that the factors of N are chosen in our definitions such that the interval [—a, +a] 
over which p(p) is non-zero is finite (i.e., of order N®) in the large N limit. 

For applications to disordered system ip is often thought of as the Hamiltonian. Its 
eigenvalues then describe the energy levels of the system. In some applications, <p is related 
to the transmission matrix/ The density of and correlation between its eigenvalues tell us 
about the conductance fluctuation in disordered metals. 

The density of eigenvalues has long been known in the literature^ to have the form 

p(n) = -P(n)y/a 2 -n 2 (1.10) 

where P(p) is a polynomial of degree 2p — 2 if the potential V is of degree 2p. The 
polynomial P(p) and the endpoint of the spectrum a depend on V in a complicated way. 
In our recent work, we have focussed on the correlation between eigenvalues. 

We have obtained the following results, > all in the large N limit. 

(1.) Using the method of orthogonal polynomials, we determine 



x <^ [cas[(p(fj)] - cas[(p(y)]][cas[N(h(ji) + (u))] + cos[N(h(p) - h(u))]] 

I (1.11) 
+ [sm[<p(p)} - sm[<p(v)]] sm[N(h(fi) + h(u))] 

t 2 

+ [sm[(p(n)\ + em[<p(v)]] sm[N(h(p) - h(v))\ ' 



The various functions /(/f), ^(a 4 ) an d h(p) are given in Ref. 5 and depend on V. Thus, 
this rather complicated expression depends on V in detail. 

(2.) With \i — v close together, of order iV -1 times some number large compared to unity 
so that there can still be a finite number of eigenvalues between p and u, and with 
both p and v at a finite distance from the end points of the spectrum, we obtain the 
universal result 

A '<"' ") ^ ^nh!^ - (L12) 

where p = ^(p — v) and 5p = \(p — v). This means that the correlation function takes 
the form 



1 



sin .c ' " 



X 



(1.13) 



x = 2rrNdpp(p). 
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(1.14) 



This result has long been known in the literature. For (p — v) small this implies that 
p(p,, v) vanishes as 

p{p^)~\^p{p)\X{p-v)\\ (1.15) 
as expected from the Van der Monde determinant in the measure. 

Note in this context that it is incorrect to replace sin 2 a; by its average 1/2, as is 
sometimes done in the literature. 

(3.) The wild oscillation of p c (p>, v) is entirely as expected since between p and v finitely 
separated there are in general O(N) eigenvalues. Thus, it is natural to smooth p c (p, v) 
by integrating over intervals dp and 8v large compared to O(N^) but small compared 
to 0(N°) centered around p and v respectively. We then obtain 5 



smooth/ n - - t^_i_ n iq) 

We find this result rather remarkable since the density p(p) is completely non-universal. 
The only dependence on V appears through a. 

(4.) We have also computed the three- and four-point connected correlation functions. 
When the oscillations in these functions are smoothed over we found that they vanish 
identically to 0(iV -3 ) and 0(iV~ 4 ) respectively. We conjectured that the smoothed 
p-point connected correlation function similarly vanishes to 0(N~ P ). 

(5.) We can show that the results in (1),(2),(3) hold for an ensemble much more general 
than the one defined in (1.1) and (1.2), namely, an ensemble defined with 

H{ V ) = ^TrV(^) + -l^TrW 1 (ip)TrW 2 (i P ) 

^ N (1.17) 

+ j^TrX 1 (i P )TrX 2 (ip)TrX 3 (i P ) + ... 

with W, X, Y, . . . etc. arbitrary polynomials. Indeed, it is easy to see 6 that this ensemble 
can be further generalized by replacing, for example, the third term in (1.13) by 

^Y, TrX iM TrX *M TrX aM (1-18) 

a 

with Xf a polynomial. This ensemble appears to us to be the most general ensemble 
invariant under unitary transformations, i.e., 

P{tf<pU) = P((p) (1.19) 
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(a 2 — pv) 



except for some rather singular examples. 

(6.) Wigner 1 also studied non- unitary ensembles, for example an ensemble of matrices 
ip whose matrix elements take on the value ±v/y/N with equal probability. Using 
a renormalization group 9 inspired approach we can show that p(p) for this class of 
matrix is universal and equal to the p{p) for the simple Gaussian ensemble defined 
with V{(f) = Vftp 2 in (1.2), namely 

p(p) = -\ja 2 -p 2 (1.20) 

that is, Wigner's well-known semi-circle law. 

(7.) In Ref. 6, we outline an argument showing that the results stated in (1), (2), (3), and 
(4) also hold for the non-unitary ensemble mentioned in (6). 

Some of our results appear to overlap with results obtained in recent literature. 10 
In particular, our smoothed universal connected correlation (1.16) appears to have been 
discovered also by Beenakker in an interesting work, 11 although we have not yet seen a 
full derivation. 

Given our correlation functions (1.11), (1.13), and (1.16), we can proceed to determine 
the mean square fluctuation of physical quantities such as conductance. For any quantity 
A((p) defined by the trace of some function of <p, the mean square fluctuation is clearly 
given by 

var A = (a 2 ^ - (A) 2 = J dpdvp c (p, v)A{p)A{v) (1.21) 

We would like to emphasize that it is incorrect to put for p c (p, v) the smoothed correlation 
p| m00th (^, v) a ^} v y 2 given in (1.16). The singularity in p™ ooth (^, v) as p - v -> 
would produce a divergent integral. Within our discussion, there is of course no difficulty 
whatsoever, since p s J ROOtil (p, v) was derived with the explicit proviso that p — v is of 0(iV ) 
and the true p c (p, v) given in (1.11) is perfectly smooth as p — v — > 0. Calculation of the 
variance of the conductance given in the recent literature 11 appears to us to involve simply 
replacing p c by ^smooth w ithout justification. 

From the definitions for p(p) and p c (p>, v), it is easy to derive the response of p(p) 
under a change in the potential V(u), as pointed out by Beenakker 11 : 

Sp(p) = -N 2 J dvp c (p, u)SV(u) (1.22) 

Again, it would be tempting to replace 11 p c (p>, v) by pf noot ^(p,, v). However, this illegiti- 
mate procedure would lead to a divergent integral. The variation 5p(p)/5V(v) appears to 
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pose a rather tedious calculation with P and a in (1.10) depending on V in a complicated 
way. 

We took V to be an even polynomial for simplicity, so that the density of eigenvalues 
is a symmetric function between the endpoints —a and +a. It is a simple matter to shift 
the spectrum. Clearly, if we replace in V by — dl (with / the unit matrix and d the 
shift) the density of eigenvalues would be non-zero between c = —a + d and b = a + d. The 
universal correlation function is then trivially shifted to read 



SmOOth ( ) = - l -~ ■ av- ■ -nr - r--; 

V ^ 2iV2 7 r 2 (^-z/) 2 [(6- /U )( /U - c )(6-z/)(z/-c)]V2 V ; 

In section II, we will develop a diagrammatic method which will enable us to study 
"time" dependent correlation between the eigenvalues, and which when the time depen- 
dence is suppressed allows us to recover many of the results mentioned above. In section 
III we compute the current-current correlation function. In section IV we study a class 
of Hamiltonians consisting of the sum of a deterministic term and a random term. The 
correlation between the eigenvalues when the deterministic term is varied is calculated. 



1 1 [-bc + \(b + c){n + v) -H 



II. "TIME" DEPENDENT CORRELATION 



In this paper, we introduce a diagrammatic method to study the correlation between 
eigenvalues for a time-dependent ensemble of hermitean matrices with a probability dis- 
tribution defined by 



2 V dt 



(2.1) 



1 P 

P((p) = - exp - / dt Tr 

Z J-T 

We take T — > oo. Here the matrix ip(t) depends on time, or more generally, on some 
external parameter we are allowed to vary. Physically, we may apply our results obtained 
below to disordered systems in which the disorder may vary. In going from (1.1) to (2.1) we 
are moving from zero-dimensional field theory (i.e., an integral) to one-dimensional field 
theory (i.e., Euclidean quantum mechanics). In the language of string theory, we move 
from a central charge c = theory to a c = 1 theory. As a byproduct we show how some 
of our previous results mentioned in Section I may be recovered as a special case. For ease 
of presentation, we will take V((p) = ^r^ 2 to be Gaussian and indicate below how our 
results may be generalized. 
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The one-point Green's function, or more generally the propagator, 



G «M = ((j^)fe> (2-2) 



can be readily determined since due to time translation invariance it does not in fact 
depend on time. Using the usual Feynman diagram expansion, we find immediately that 
in the large N limit Gij(z) is given by planar Feynman diagrams (generalized rainbow) as 
indicated in Fig. (la). 

It is perhaps useful, borrowing the terminology of large N QCD from the particle 
physics literature, to speak of the single line in Fig (la) as representing quark propagators 
and the double lines as gluon propagators. The quark propagator is given simply by 1/z 
while the gluon propagator is given by 

DijM*) = (^y(*)v«(0)> 

1 fdu . _ 1 _ m m (2-3) 



= Wjkjj J ^^2^2 = Mjfc^e 

We will now immediately generalize to arbitrary time dependence by replacing (2.3) 

by 

2 

D l]M (t) = 6 a 6 jk ^e-«® (2.4) 

For instance, the time dependence in (2.1) may be changed in such a way that the "mo- 
mentum space" propagator u; 2 + m 2 in (2.4) is replaced by bJ ' 1 +'y\uj\+Tri i or ZF+^F+m 1 ^ or 
example. In other words, in (2.1) we can replace \{^) 2 by ipK(^)ip with K any reason- 
able function. Our only requirement is that u is a smooth function of t and does not blow 
up as t goes to zero. 

Introducing as usual the one-particle irreducible self energy Ej ? -(z), we can write the 
generalized rainbow integral equation (Fig. (lb)) 

E(z) = ^j^i) = °* G W ( 2 - 5 ) 

Here we have used the fact, immediately obvious from examining the Feynman diagrams, 
that ^ij(z) is equal to 8{j £(2:) and the fact that the gluon propagator only appears at 
equal time ^(0) = S^Sj^^a 2 . Note that the quark does not know about time. Solving 
the quadratic equation for E we obtain the Green's function as defined in (1.3) 

G(z) = ^(z-Vz 2 -^) (2.6) 



Taking the absorptive part we recover immediately Wigner's semicircle law as given in 
(1.20). 

Incidentally, within this diagrammatic approach, band matrices can be treated im- 
mediately. Let the matrices (p be restricted so that (pij vanishes unless \i — j\ < bN/2 
with b < 1. Such matrices describe, for example, the hopping of a single electron on a 
one-dimensional lattice with random hopping amplitudes. The essential feature is that 
from each site the electron can hop to 0(1/N) sites. Looking at the Feynman diagrams, 
we see that in the generalized rainbow integral equation we simply restrict the range of 
summation from N to bN and thus instead of (2.3) we obtain 

E(z) = ba 2 —^ = ba 2 G(z) (2.7) 

Thus, we have the same distribution of eigenvalues with a suitable re-definition of the 
endpoints. 

Let us now move on to the connected 2-point Green's function 

/ 1 1 1 

G c (z,w,t) = ( — T,tr —tr 



N 2 z-(p(t) w-(p(0)/ c 

1 oo oo (2-8) 

m=0 n=0 

Henceforth for the sake of notational clarity we will set a to unity; it can always be 
recovered by dimensional analysis. Diagrammatically, the expression for G c (z, w, t) can 
be described as two separate quark loops, carrying "momentum" z and w respectively, 
interacting by emitting and absorbing gluons (see Fig. (2a)). 

With a Gaussian distribution for (p, we can readily "Wick-contract" the expression 
{trip m (t)tnp n (0))(j. Let us begin by ignoring contractions within the same trace (in which 
case m and n are required to be equal). In the large N limit, the dominant graphs (see 
Fig. (2b)) are given essentially by "ladder graphs" (with one crossing) which immediately 
sum to 

N 2 G c {z,w) = J^^J^—J-^ (2.9) 

V ' \ ZW> 

We next include Wick-contractions within the same trace in (tr(p m tr(p n ). We see that 
graphically these contractions describe vortex and self energy corrections. The vortex 
corrections can be immediately summed: the expression in (2.9) is to be multiplied by two 
factors, the factor 

(2.10) 
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and a similar factor with z replaced by w. Finally, according to our earlier discussion, 
self-energy corrections are included immediately by replacing the bare quark propagator 
1/z by the dressed propagator G(z) (and similarly for 1/w of course.) We obtain finally 
the remarkably compact result 



N 2 G c (z,w) 



1 



(1 - G(z)G(w)y 



1-G 2 (z)_ 



G 2 {w) 



1-G 2 ( 



w 



(2.11) 



Finally, we have to put in the time dependence, which we have ignored so far. We note 
first of all that the vertex and self-energy corrections contain no time dependence, since 
there the gluons always begin and end on the same quark line. Thus, the time dependent 
2-point connected Green's function can be written down immediately as 



N 2 G c (z,w,t) 



-u{t) 



(l-e-< t )G{z)G{w)) 2 



1 - G 2 {z) 



G 2 (w) 
1-G 2 (w) 



(2.12) 



To obtain the connected correlation function between energy eigenvalues, we have to 
take the double absorptive part of G c (z,w,t) as indicated in (1.7). It is most convenient 
to introduce angular variables: from (2.6) we see that we may write 



where rj = the sign of e and 

Similarly, we write 

with £ = the sign of 5 and 



G(/i + ie) = -irje^ 6 

sin 9 = \xja. 
G(v + iS) = -iZe** 



(2.13) 

(2.14) 
(2.15) 
(2.16) 



sin</> = v la. 

As /i and v vary over their allowed ranges, from —a to +a, 9 and (f> range vary from — tv/2 

to 7T/2. 

We can now readily compute (in the notation of (1.8)) 



8N 2 G C (++) = 



1 



cos 9 cos 4> I [1 + ch u cos(6' + (p) — ish u sin(# + 
Proceeding in this way, we obtain one of our main results 

- 16n 2 N 2 p c (v, u) 

( 1 + ch ucos(6> + 4>) 



(2.17) 



cos#cos(/> \ [ch u + cos(# + 4>)] 2 
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+ ((f) -> -(f) + 7r) 



(2.18) 



Note that crossing symmetry [ji <-> v, t — > — £) clearly holds. For £ = 0, that is, n = 
(since any non-zero it(0) can be absorbed into cr 2 ), we recover immediately our previous 
result (1.16), obtained by the orthogonal polynomial method. Note that time acts as a 
regulator for the singularity when fj, = v: for time not equal to zero, we can set \i = v 
without difficulty and obtain 

8n 2 N 2 p c (ia, //) = — + • • • (2-19) 
cos z u l 



For ji 7^ v and small time u « 6 — </> we have 

1 _ 1 ms ffl 4- rh) „ 

+ (</> 



-167r 2 iV 2 p c (At, i/) ' 1 



cos 6* cos 11 + cos(6> + 4>) 



l l-^cos(fl + 0) ^ 2 
1 + cos(6> + 0) 




In the long time limit, w — > oo, we find 

47r 2 A^ 2 p c (^, i/) -> e _u tan^tan0 (2.21) 

As expected the correlation vanishes exponentially in time. Notice however that a 
memory of the spatial correlation is retained even at arbitrarily large time. 

We thus conclude that the density-density correlation is universal in space for all time, 
in the sense that it does not depend on V at all. 

It is noteworthy that in the diagrammatic approach the correlation function is "auto- 
matically" smoothed. An interchange of limits is responsible. Here in computing G(z, w) 
we are taking to infinity and then letting z approach the real axis. In our previous 
work 5 we use orthogonal polynomials to calculate p c (fJ>, v) directly. In effect, we let z sit 
on the real axis before taking N to infinity and thus the discrete pole structure is visible. 
In this sense, the orthogonal polynomial method is more powerful and informative. 

Having now analyzed the Gaussian case we now discuss how our results could hold 
more generally. We distinguish between the "trace class" defined in (1.1) and the "Wigner 
class" defined in (6). 

For the Wigner class, let us focus on the example in which the probability of the 
distribution matrix element Lpij is given by 

P(<Pij) oc e _Ar2(l ^ |2 "^ )2 (2.22) 

corresponding essentially to the example mentioned in item (6) in Section I (see Ref. 6). 
The quartic interaction ~ A 2 |^^ | 4 would contribute to Feynman diagrams such as the 
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one in fig (3a). Counting powers, we see that this graph is of order A~ 4 AA 2 = A -1 
and so is suppressed relative to the graphs in Fig. (la). Reasoning along this line, we see 
immediately that the distribution of eigenvalues is universal, a long-known result that we 
also derived recently using a renormalization group inspired approach. 6 As a bonus, we 
obtain immediately in the present diagrammatic approach that the correlation function is 
also universal. (Incidentally, this result is not at all easy to obtain with the renormalization 
group approach of Ref. 6.) 

For the trace class, the interaction, for example the quartic term in V((p), would 
generate Feynman diagrams such as the one in figure (3b). Counting powers of A we see 
that this graph is in no way suppressed relative to the graphs in Fig. (la). This is in fact 
a gratifying conclusion as we know from Ref. 8 (See eq. (1.10)) that in the trace class, in 
sharp contrast to the Wigner class, the distribution of eigenvalues is in fact not universal. 
It appears to us that within the diagrammatic approach, it would be rather involved to 
demonstrate the universality of the correlation function for the trace class, namely the 
result we obtained in Ref. 5 using the method of orthogonal polynomials. We would have 
to show that the effects of the arbitrary polynomial interactions contained in V in (1.1) 
can be summed up and absorbed completely in the endpoint value a of the spectrum. 

We find it remarkable that in this subject results easily obtained in one approach are 
apparently rather difficult to prove in another. 



As is well-known, the ensemble in (1.1) may be thought of as describing the statistical 
mechanics of a gas in 1 + dimensional space time. The partition function 



where C is an irrelevant overall constant. The logarithmic term comes from the well-known 
Jacobian connecting dip to d\\d\2 ■ ■ .dXpf. Regarding as the position of "z th particle" 
on the real line, we see that (3.1) describes a one-dimensional gas of A particles interacting 
with each other via a logarithmic repulsion while confined by an external potential V(X). 
Intuitively then, it becomes entirely clear that the density of the gas p(fi) has no reason 



III. CURRENT-CURRENT CORRELATIONS 




(3.1) 
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to be universal: it should certainly depend on V. It is less clear why the change in density 
Sp(p) at p due to a change in the potential 8V(v) at v should be universal, and indeed, 
this universality holds only when we smooth over the discrete character of the gas. 

The generalization in (2.1) then corresponds to allowing the particles to move in a 1 + 1 
dimensional spacetime. With the density operator defined by 

PGM) = ^J>Gu-A;(t)) 



we clearly have the conservation law 



with the current operator defined by 



dt dp, 



J(M) = ^£^(^-A 2 (t)) 

i 

Thus, we are led to study the current- current correlation function 
Note that (J J) is connected by definition. 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



The double integral in (3.5) may be explicitly evaluated. We find that in effect p c (p, v) 
may be written as 



-4n 2 N 2 p c (p, u) : 
In particular, at equal time, we have 



d_d_ 
dp dv 



log 



ch u + cos(6> + <fi) 
ch u — cos(6> — 4>) 



Ait 2 N 2 Pc (p, v) = log 



d d 



dp dv 



a 2 - pv + A/(a 2 - p' 2 )(a' 2 - v 2 ) 
a 2 — pv — a/ (a 2 — p 2 )(a 2 — v 2 ) 



(3.6) 



(3.7) 



An equivalent form reads 



2iv 2 N 2 p c (p, v) = ^tMo 
dp dv 



d d 



a—fj, 
a+/i 



a+is 



a— M I / a— v 
a+fi ' Y a+v _ 



(3.7a) 



The current- current correlation function then follows immediately 



4n 2 N 2 (J(p, t)J(v, 0)) = ^log 



d 2 



dt 2 



ch u + cos(6> + 4>) 
ch u — cos((6> — </>) 



sh u 



-u + 



(1 + ch u cos(6> + 



ch u + cos(8 + 4>) (ch u + cos(# + 4>)) 2 
12 



U — ((f) — > —(f) + 7r) 

(3.8) 



We have thus obtained the current-current correlation function for arbitrary separation 
in space and time. 

First, the dependence on space is universal: as a function of 6 and 0, the current- 
current correlation, just like the density- density correlation from which it is derived, does 
not depend on the potential V. 

The dependence on time, in contrast, is no n- universal: clearly, the dependence of u on 
t enters. In special cases, however, the specific functional form may be seen to drop out. 
From (3.8), we see that at the same point in space, that is when 9 = (f), and u small, we 
have 

8 2 

2n 2 N 2 (J(p, t)J(v, 0)> = ^ log u + . . . (3.9) 

Thus, if as t — > 0, u vanishes like u — > at 7 , then 

2n 2 N 2 (JJ) -> - 7 /t 2 (3.10) 

We have universality in the sense that the unknown constant a has dropped out. With the 
further assumption that 7 = 1, which is reasonable but certainly not required, we obtain 
the universal statement 

2n 2 N 2 (JJ) -> -1/t 2 (3.11) 

We also observe the curiosity that at 9 = = ±7r/2, the current-current correlation 
vanishes identically for all time. 

For large time, we obtain 

7i 2 N 2 (J(/j, t)J(u, 0)) -> cosflcos^-u 2 - u)e~ u (3.12) 

This is of course not independent of (2.19). 

We have learned that the universality of (J J) at small time has already been derived 
by Szafer and Altschuler 12 and by Beenakker 13 using apparently rather different methods 
and implicitly assuming that 7 = 1. To our knowledge, the complete form of {J J) in (3.8) 
has not appeared before in the literature. 

We can immediately generalize the preceding discussion to the case where many ex- 
ternal parameters t\, t<i-, ■ ■ ■ , . . . are varied. This may be described picturesquely as a 
many-time world in which the conservation laws 

dJ k dp 



13 



(3.13) 



hold where the current with respect to time is defined by 

f)\ ■ 

J k (ri = J EdT 5{fi - Xi{t)) (3 - 14) 

i k 

The relation (3.5) immediately generalizes to 

(J k (p, t)J l {^ 0)> = -d k d l f dp' f dv'(p(p', t)p{v\ 0)> c (3.15) 

J — a J —a 

where d k = For the special case where u is a function of t = (^ fc ) 1//2 j f° r example, 
we obtain easily that 

2^N*(j k (p, t)J l {^ 0)) - (3.16) 
for small time, with an assumption similar to the one that leads to (3.11). 



IV. DETERMINISTIC PLUS RANDOM 



Our diagrammatic approach allows us to study immediately the eigenvalues of a Hamil- 
tonian of the form H — Hq + <p which consists of the sum of a deterministic piece Hq and 
a random piece <p with a probability distribution such as in (1.1). Pastur 14 has found the 
interesting relation that 

G(z) = G (z - G(z)) (4.1) 

where G and Gq are the Green's functions for H and Hq respectively. We now show that 
Pastur's relation follows immediately from our diagrammatic analysis. 

Let us consider the Gaussian case and let Hq be diagonal with diagonal elements q. 
Looking at the relevant Feynman diagrams, we see that we simply have to replace the 
inverse quark propagator z by z — ej and thus obtain 

We recognize this as (4.1). The relation (4.1) while interestingly compact is not terribly 
useful in practice as in solving for G(z) one would encounter a polynomial equation of 
degree N + 1. 

14 



We will now demonstrate the power of the diagrammatic approach by showing how 
we can immediately go beyond Pastur's relation and study correlation. Consider the 
following class of physical problems. Suppose we change some external parameters so 
that the deterministic Hamiltonian Hq is changed to Hq. We would like to compute the 
correlation between the spectra of H and H 1 = Hq + (p, in other words, we would like to 
compute 

G c {z,w,Hq,H'q) = (hr i hr * \ (4.3) 

\N z- Hq- <pN w- H Q -<p/ c 

where the average, as before, is over the distribution of (p. 

An example of this class of problems was recently studied by Simons and Altschuler 15 . 
They considered the problem of a single non-interacting electron moving in a ring threaded 
by a magnetic flux and with the electron scattering on impurities in the ring. The magnetic 
flux is then changed to some other value with H changed accordingly to H' . The correlation 
between the spectra of H and H' is apparently of great interest in the physics of mesoscopic 
systems. 

We see immediately from Fig. (2b) that using the diagrammatic approach we can 
determine G c (z,w, Hq, Hq) quite readily. Let us define 

where is the solution of (4.2). We assume that Hq is also diagonal, with diagonal 
elements e'. We define the analog of gi(z) for Hq, namely 

hi(w) = , 1 (4.5) 

w — — Ij'(w) 

where T,'(w) is the solution of the analog of (4.2) for Hq, namely (4.2) with q — > e[, z — > w, 
and E(z) — > Y,'(w). Let us also introduce the short hand notation 

9-9= ^i9 % 

g-h = —X igi hi (4.6) 



9 ■ h = jf^i9i h i 



and so forth. Then 



G c(^H ,H Q ) = ( { l- g . h f \-g-gl-h-h (4 ' 7) 
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We see that this collapses to (2.11) when Hq = Hq = 0. 

The number of physical situations covered by the result of this section is very large. In 
particular, it should allow us to verify the intuitive expectations concerning the universality 
of the correlations: the general belief is that if the energy scale is such that the system 
can explore the full extent of the disorder, one should recover the correlations for the pure 
random matrix correlations, irrespective of the non-random Hq. It is far from obvious 
from the explicit representation (4.7). Clearly, our result (4.7) could also be used to study 
questions such as localizations, or the influence of white noise in various physical situations. 
Thus, we believe that the result in (4.7) would prove to be of importance in studying many 
disordered systems. 

APPENDIX 

We have focussed on random hermitean matrices. In some physical situations, the 
Hamiltonian matrix is in fact real symmetric. It is well-known in the literature (see for 
example Ref. 7) that going from the case of hermitean random matrices to the case of real 
symmetric matrices we simply insert in various formulas appropriate factors of 2. This 
is most easily seen by looking at (3.1): for (p real symmetric matrices we would have a 
factor of 1 /2 in front of the logarithmic repulsion. This comes about because the Jacobian 
connecting d<p to d\id\2 . . .dXjy for the hermitean case is the positive square root of the 
corresponding Jacobian for the real symmetric case. 

Let us now sketch exceedingly briefly how this factor of 2 emerges in the diagrammatic 
approach. Typically, we encounter (trip n trip n ) . To indicate how the argument goes, let 
us consider only Wick-contractions between the two traces. After the first contraction, we 
have 

n{vijVaf3){{v n ~ l )ji{v n ~ l )p a ) = «(* w ^+^* ja )((^ 1 ) J! (^ 1 )/J a ) = 2nS 6 ja ^r 1 

(Al) 

The last line follows from the fact that ip is symmetric. Proceeding, we find easily that the 
above is equal to 2nN n . Note that for hermitean matrices the expression in the parenthesis 
would read Si a Sj/3 instead, and this accounts for the factor of 2 alluded to above. 

Graphically, the fact that (<Pij<p a p) is equal to Si a djp+d^dj a rather than Si a 5jp means 
that the double lines in the gluon propagator in the Feynman diagrams can be twisted. 
The direct of counting of graphs becomes considerably more involved. 
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We also note that the method of orthogonal polynomials used in Ref. 5 becomes rather 
complicated when we deal with real symmetric matrices. 
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FIGURE CAPTIONS 

Fig (l.a) Feynman diagram expansion for the Green's function G(z). 

Fig (l.b) The generalized rainbow equation for the one-particle-irreducible self energy 
E(z). 

Fig (2. a) Feynman diagram expansion for the connected two-point Green's function. 
The dotted lines here represent the "gluon" propagator. 

Fig (2.b) Some typical low order Feynman diagrams contributing to the connected 
two-point Green's function. 

Fig (3) Non-Gaussian corrections to the propagator in the (a) Wigner class and in the 
(b) trace class. 
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